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Abstract 


BACKGROUND: The future of the SARS-CoV-2 pandemic hinges on virus evolution and 
duration of immune protection of natural infection against reinfection. We investigated duration 
of protection afforded by natural infection, the effect of viral immune evasion on duration of 
protection, and protection against severe reinfection, in Qatar, between February 28, 2020 and 


June 5, 2022. 


METHODS: Three national, matched, retrospective cohort studies were conducted to compare 
incidence of SARS-CoV-2 infection and COVID-19 severity among unvaccinated persons with a 
documented SARS-CoV-2 primary infection, to incidence among those infection-naive and 


unvaccinated. Associations were estimated using Cox proportional-hazard regression models. 


RESULTS: Effectiveness of pre-Omicron primary infection against pre-Omicron reinfection was 
85.5% (95% CI: 84.8-86.2%). Effectiveness peaked at 90.5% (95% CI: 88.4-92.3%) in the 7" 
month after the primary infection, but waned to ~70% by the 16" month. Extrapolating this 
waning trend using a Gompertz curve suggested an effectiveness of 50% in the 22™ month and 
<10% by the 32™ month. Effectiveness of pre-Omicron primary infection against Omicron 
reinfection was 38.1% (95% CI: 36.3-39.8%) and declined with time since primary infection. A 
Gompertz curve suggested an effectiveness of <10% by the 15" month. Effectiveness of primary 
infection against severe, critical, or fatal COVID-19 reinfection was 97.3% (95% CI: 94.9- 
98.6%), irrespective of the variant of primary infection or reinfection, and with no evidence for 


waning. Similar results were found in sub-group analyses for those >50 years of age. 


CONCLUSIONS: Protection of natural infection against reinfection wanes and may diminish 


within a few years. Viral immune evasion accelerates this waning. Protection against severe 
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reinfection remains very strong, with no evidence for waning, irrespective of variant, for over 14 


months after primary infection. 
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Introduction 


The future of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic is 
uncertain, but it hinges on virus evolution and the level and duration of immune protection of 
natural infection against reinfection. While current coronavirus disease 2019 (COVID-19) 
vaccines had a critical role in reducing COVID-19 hospitalizations and deaths, their rapidly 
waning immune protection, particularly against the Omicron (B.1.1.529) variant,** limits their 
role in shaping the future of SARS-CoV-2 epidemiology compared to other vaccines, such as 


vaccinia, which eradicated smallpox.’ 


Seasonal “common-cold” coronaviruses are characterized by short-term immunity against mild 
reinfection,'® but long-term immunity against severe reinfection.” While SARS-CoV-2 infection 
with the original virus or pre-Omicron variants elicited >80% protection against reinfection with 
the original virus!!? or with Alpha (B.1.1.7),!* Beta (B.1.351),!° and Delta (B.1.617.2)!® 
variants, protection against reinfection with Omicron subvariants is below 60%.'°!” Reinfections 


have become common since Omicron emergence." 


We sought to answer three questions of relevance to the future of this pandemic: 1) When 
infected with a pre-Omicron variant, how long does protection persist against reinfection with 
pre-Omicron variants? 2) When infected with a pre-Omicron variant, how long does protection 
persist against reinfection with an Omicron subvariant? 3) When infected with any variant, how 
long does protection persist against severe, critical, or fatal COVID-19? Answers to these 
questions help us to understand duration of protection resulting from natural-infection, effects of 
viral evasion of the immune system on this duration, and effectiveness of natural infection 


against COVID-19 severity when reinfection occurs. 
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Three studies were conducted to answer these questions in Qatar, a country that experienced five 
SARS-CoV-2 waves dominated by each of the original virus,!! Alpha,'* Beta,!> Omicron BA. 1 
and BA.2,'* and currently Omicron BA.4 and BA.5, in addition to a prolonged low-incidence 


phase dominated by Delta.* 
Methods 
Study population and data sources 


This study was conducted in the population of Qatar and analyzed COVID-19 data for laboratory 
testing, vaccination, hospitalization, and death, retrieved from the national digital-health 
information platform. Databases include all SARS-CoV-2-related data, with no missing 
information since pandemic onset, such as all polymerase chain reaction (PCR) tests, and starting 
from January 5, 2022, rapid antigen tests conducted at healthcare facilities. Further descriptions 


of the study population and these national databases were reported previously. *15-17-19-20 


Study designs and cohorts 


We conducted three matched, retrospective cohort studies that emulated randomized “target” 
trials.?°?! In each study, incidence of infection or of severe,” critical,” or fatal? COVID-19 was 
compared in the national cohort of individuals with a documented SARS-CoV-2 primary (first) 
infection prior to vaccination (designated the primary-infection cohort) to the national (control) 
cohort of individuals who are infection-naive and unvaccinated (designated the infection-naive 


cohort). 


Documentation of infection in all cohorts was based on positive PCR or rapid antigen tests. 


Laboratory methods are found in Supplementary Appendix Section S1. Classification of 
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COVID-19 case severity (acute-care hospitalizations),”” criticality (intensive-care-unit 


2 
ed 


hospitalizations),”” and fatality” followed World Health Organization guidelines (Section S2). 


Cohort matching and follow-up 


Individuals in the primary-infection cohort were exact-matched in a one-to-one ratio by sex, 10- 
year age group, nationality, and comorbidity count (none, 1-2 comorbidities, 3 or more 
comorbidities) to individuals in the infection-naïve cohort, to control for differences in risk of 
SARS-CoV-2 infection in Qatar. !?4+? Matching by these factors was shown to provide adequate 
control of differences in risk of infection.47*3! Matching was also done by the calendar week of 
SARS-CoV-2 testing. That is, an individual who was diagnosed with a primary infection in a 
specific calendar week was matched to an infection-naive individual who had a record of a 
SARS-CoV-2-negative test in that same week (Figures $1-S3). This matching ensures that all 
individuals in all cohorts had active presence in Qatar at the same calendar time. Matching was 
performed using an iterative process so that each individual in the infection-naive cohort was 
alive, infection-free, and unvaccinated at the start of follow-up. 

SARS-CoV-2 reinfection is conventionally defined as a documented infection >90 days after an 
earlier infection, to avoid misclassification of prolonged PCR positivity as reinfection. '7'*-'© 
Therefore, each matched pair was followed from the calendar day an individual in the primary- 


infection cohort completed 90 days after a documented primary infection. 


For exchangeability, both members of each matched pair were censored on the date of first-dose 
vaccination of an individual in either cohort.?°? Individuals were followed up until the first of 
any of the following events: a documented SARS-CoV-2 infection, i.e., the first PCR-positive or 
rapid-antigen-positive test after the start of follow-up, regardless of symptoms, or first-dose 


vaccination (with matched pair censoring), or death, or end of study censoring. 
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Pre-Omicron Reinfection Study 


This study estimated effectiveness of a pre-Omicron primary infection against reinfection with a 
pre-Omicron variant. Any individual with a documented primary infection between February 28, 
2020 (pandemic onset in Qatar) and November 30, 2021 was eligible for inclusion in the 
primary-infection cohort, provided that the individual received no vaccination before the start of 
follow-up, 90 days after primary infection. Any individual with a SARS-CoV-2-negative test 
during this period was eligible for inclusion in the infection-naive cohort, provided that the 
individual had no record of infection or vaccination before the start of follow-up. Follow-up was 
from the 90" day after primary infection until November 30, 2021 (first evidence of Omicron in 
Qatar!®° to ensure that incidence during the study was only due to a pre-Omicron variant). The 
primary study outcome was incidence of infection. The secondary outcome was incidence of 


severe, critical, or fatal COVID-19. 
Omicron Reinfection Study 


This study estimated effectiveness of a pre-Omicron primary infection against reinfection with 
an Omicron subvariant. Any individual with a documented primary infection from February 28, 
2020 until November 30, 2021 was eligible for inclusion in the primary-infection cohort, absent 
any record of reinfection or vaccination before the start of follow-up. Any individual with a 
SARS-CoV-2-negative test during this period was eligible for inclusion in the infection-naive 
cohort, absent any record of infection or vaccination before the start of follow-up. Follow-up was 


from December 19, 2021 (onset of the Omicron wave in Qatar), 16? 


if the primary infection 
occurred >90 days before this date. Follow-up was from the 90" day after primary infection if 


the primary infection occurred <90 days before December 19, 2021. Follow-up was until June 5, 
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2022. The primary study outcome was incidence of infection. The secondary outcome was 


incidence of severe, critical, or fatal COVID-19. 


COVID-19 Severity Reinfection Study 


This study estimated effectiveness of a primary infection against severe, critical, or fatal 
COVID-19 reinfection, irrespective of the variant of primary infection or reinfection. Any 
individual with a documented primary infection between February 28, 2020 and June 5, 2022 
was eligible for inclusion in the primary-infection cohort, provided that the individual received 
no vaccination before the start of follow-up, 90 days after primary infection. Any individual with 
a SARS-CoV-2-negative test during this period was eligible for inclusion in the infection-naive 
cohort, absent any record of infection or vaccination before the start of follow-up. The primary 
study outcome was incidence of severe, critical, or fatal COVID-19. The secondary outcome was 


incidence of infection. 


Statistical analysis 


Eligible and matched cohorts were described using frequency distributions and measures of 
central tendency, and were compared using standardized mean differences (SMDs). An SMD 
<0.1 indicated adequate matching. Cumulative incidence of infection (defined as the proportion 
of individuals at risk, whose primary endpoint during follow-up was a reinfection for the 
primary-infection cohort, or an infection for the infection-naive cohort) was estimated using the 
Kaplan—Meier estimator method. Incidence rate of infection in each cohort, defined as the 
number of identified infections divided by the number of person-weeks contributed by all 
individuals in the cohort, was estimated with its 95% confidence interval (CI) using a Poisson 


log-likelihood regression model with the Stata 17.0 stptime command. 
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The hazard ratio, comparing incidence of infection in both cohorts, and the corresponding 95% 
CI were calculated using Cox regression adjusted for matching factors with the Stata 17.0 stcox 
command. Schoenfeld residuals and log-log plots for survival curves were used to test the 
proportional-hazards assumption and to investigate its adequacy. 95% CIs were not adjusted for 
multiplicity; thus, they should not be used to infer definitive differences between cohorts. 
Interactions were not considered. Effectiveness against reinfection was estimated using the 


equation: Effectiveness = 1 — adjusted hazard ratio. Analogous analyses were used when the outcome 


was severe, critical, or fatal COVID-19. 


Subgroup analyses were conducted to investigate waning of protection over time. Adjusted 
hazard ratios were estimated by month since primary infection using separate Cox regressions 
with "failures" restricted to specific months, in the Pre-Omicron Reinfection Study and the 
COVID-19 Severity Reinfection Study. Adjusted hazard ratios were also calculated in the 
Omicron Reinfection Study, but stratified by 3-calendar-month primary-infection sub-cohorts. 
Additional analyses restricting matched cohorts to those >50 years of age were conducted. 
Sensitivity analyses adjusting effectiveness estimates for differences in testing frequency 


between cohorts were conducted. 


Waning of protection was fitted to the Gompertz function* using the Stata 17.0 nl command. 
This function has been used to describe decay of immunity, such as against smallpox,’ and 
provided a suitable description of waning of protection as assessed by an empiric goodness-of- 
fit. Statistical analyses were conducted using Stata/SE version 17.0 (Stata Corporation, College 


Station, TX, USA). 


Oversight 
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Hamad Medical Corporation and Weill Cornell Medicine-Qatar Institutional Review Boards 
approved this retrospective study with a waiver of informed consent. The study was reported 
following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) 


guidelines. The STROBE checklist is found in Table S1. 
Results 
Pre-Omicron Reinfection Study 


Figure S1 shows the population selection process and Table 1 describes baseline characteristics 
of the full and matched cohorts. The matched cohorts each included 290,638 individuals. The 
study was conducted on the total population of Qatar, and thus the study population is 
representative of the internationally diverse, but predominantly young and male population of 


Qatar (Table S2). 


There were 1,806 reinfections in the primary-infection cohort during follow-up, of which 6 
progressed to severe and 1 to fatal COVID-19 (Figure S1). Meanwhile, there were 11,957 
infections in the infection-naive cohort, of which 297 progressed to severe, 19 to critical, and 12 
to fatal COVID-19. Cumulative incidence of infection was 1.7% (95% CI: 1.6-1.8%) for the 
primary-infection cohort and 9.6% (95% CI: 9.4-9.9%) for the infection-naive cohort, 15 months 


after the primary infection (Figure 1A). 


The overall hazard ratio for infection, adjusted for sex, 10-year age group, 10 nationality groups, 
comorbidity count, and SARS-CoV-?2 test calendar week, was 0.14 (95% CI: 0.14-0.15; Table 3). 
Effectiveness of pre-Omicron primary infection against pre-Omicron reinfection was 85.5% 
(95% CI: 84.8-86.2%). Effectiveness increased slowly after the primary infection and reached 


90.5% (95% CI: 88.4-92.3%) in the 7" month after the primary infection (Figure 2A). Starting in 


10 
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the 8" month, effectiveness waned slowly and reached ~70% by the 16" month. Fitting the 
waning of protection to a Gompertz curve suggested that effectiveness reaches 50% in the 224 


month and <10% by the 32" month (Figure 3). 


Effectiveness of pre-Omicron primary infection against severe, critical, or fatal COVID-19 due 
to pre-Omicron reinfection was 98.0% (95% CI: 95.7-99.0%; Table 3). In the additional analysis 
restricting the matched cohorts to the sub-cohorts >50 years of age (25,595 individuals), 
effectiveness against reinfection and against severe, critical, or fatal COVID-19 reinfection was 
90.7% (95% CI: 88.4-92.5%) and 97.4% (95% CI: 91.9-99.2%), respectively. In the sensitivity 
analysis adjusting the overall hazard ratio by the ratio of testing frequency, effectiveness against 


reinfection was 79.5% (95% CI: 78.4-80.5%). More results are in Section S3. 
Omicron Reinfection Study 


Figure S2 shows the process of population selection and Table 1 describes the baseline 
characteristics of the full and matched cohorts. The matched cohorts each included 120,483 


individuals. The cohorts are representative of Qatar’s population (Table S2). 


There were 7,995 reinfections in the primary-infection cohort during follow-up, of which 5 
progressed to severe COVID-19 (Figure S2). Meanwhile, there were 12,230 infections in the 
infection-naive cohort, of which 26 progressed to severe, 7 to critical, and 5 to fatal COVID-19. 
Cumulative incidence of infection was 6.8% (95% CI: 6.7-6.9%) for the primary-infection cohort 
and 10.4% (95% CI: 10.2-10.6%) for the infection-naive cohort, 165 days after the start of 


follow-up (Figure 1B). 


The overall adjusted hazard ratio for infection was 0.62 (95% CI: 0.60-0.64; Table 3). 


Effectiveness of pre-Omicron primary infection against Omicron reinfection was 38.1% (95% 


11 
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CI: 36.3-39.8%). Effectiveness varied for the primary-infection sub-cohorts (Figure 2B). It was 
~60% for those with a more recent primary infection, between June 1, 2021 and November 30, 


2021, during Delta-dominated incidence.**+*> 


Effectiveness declined with time since primary 
infection and was 17.0% (95% CI: 10.1-23.5%) for those with a primary infection between 
December 1, 2020 and February 28, 2021, during Alpha-dominated incidence.*3435 However, 
higher effectiveness of ~50% was estimated for those with a primary infection before August 31, 
2020, during original-virus incidence (note discussion in Section $3).+*4*> Fitting the waning of 


protection to a Gompertz curve suggested that effectiveness reaches 50% in the 8" month after 


primary infection and <10% by the 15" month (Figure 3). 


Effectiveness of pre-Omicron primary infection against severe, critical, or fatal COVID-19 due 
to Omicron reinfection was 88.6% (95% CI: 70.9-95.5%; Table 3). In the additional analysis 
restricting the matched cohorts to the sub-cohorts >50 years of age (6,304 individuals), 
effectiveness against reinfection and against severe, critical, or fatal COVID-19 reinfection was 
21.6% (95% CI: 11.1-31.0%) and 84.6% (95% CI: 59.7-94.1%), respectively. In the sensitivity 
analysis adjusting the overall hazard ratio by the ratio of testing frequency, effectiveness against 


reinfection was 31.7% (95% CI: 29.7-33.6%). More results are in Section S3. 
COVID-19 Severity Reinfection Study 


Figure S3 shows the process of population selection and Table 2 describes baseline 
characteristics of the full and matched cohorts. The matched cohorts each included 407,214 


individuals, representative of Qatar’s population (Table S2). 


There were 7,082 reinfections in the primary-infection cohort during follow-up, of which 9 
progressed to severe and 1 progressed to fatal COVID-19 (Figure S3). Meanwhile, there were 


21,645 infections in the infection-naive cohort, of which 315 progressed to severe, 25 to critical, 
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and 18 to fatal COVID-19. Cumulative incidence of severe, critical, or fatal COVID-19 was 
0.00% (95% CI: 0.00-0.01%) for the primary-infection cohort and 0.21% (95% CI: 0.19-0.23%) 


for the infection-naive cohort, 15 months after the primary infection (Figure 1C). 


The overall adjusted hazard ratio for severe, critical, or fatal COVID-19 was estimated at 0.03 
(95% CI: 0.01-0.05; Table 3). Effectiveness of primary infection with any variant against severe, 
critical, or fatal COVID-19 due to reinfection with any variant was 97.3% (95% CI: 94.9- 
98.6%). Variation by month after primary infection was negligible, with no evidence for waning 


(Figure 3C). Effectiveness was ~100% up to the 14" month after primary infection. 


Effectiveness of primary infection with any variant against reinfection with any variant was 
69.4% (95% CI: 68.6-70.3%; Table 3). In the additional analysis restricting the matched cohorts 
to the sub-cohorts >50 years of age (31,108 individuals), effectiveness against reinfection and 
against severe, critical, or fatal COVID-19 was 75.3% (95% CI: 72.0-78.2%) and 95.4% (95% 


CI: 89.4-98.0%), respectively. More results are in Section S3. 
Discussion 


Protection of natural infection waned with time after primary infection, prior to Omicron 
emergence, and reached ~70% by the 16" month. This waning likely reflects genuine waning in 
biological immunity rather than viral immune evasion, as pre-Omicron variants demonstrated 
much less immune evasion than Omicron.!*'° This waning in natural immunity mirrors that of 
vaccine immunity,*°*° but at a slower rate. Vaccine immunity may last for only a year,*°*° but 
natural immunity, assuming Gompertz decay, may last for 3 years, as also suggested by long- 
term follow-up of SARS-CoV-1-associated antibodies,*° and incidentally not dissimilar to 


pandemic-influenza-associated antibodies.” 
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Immune evasion of Omicron subvariants reduced overall protection of pre-Omicron natural 
immunity and accelerated its waning (Figure 3), mirroring the effect of Omicron on vaccine 
immunity, but at a slower rate. Vaccine immunity against Omicron subvariants lasts for <6 
months,*”* but pre-Omicron natural immunity, assuming Gompertz decay, may last for just over 


a year. 


Despite waning protection against reinfection, strikingly, there was no evidence for waning of 
protection against severe COVID-19 at reinfection. This remained ~100%, even 14 months after 
the primary infection, with no appreciable effect for Omicron immune evasion in reducing it. 
This pattern also mirrors that of vaccine immunity, which wanes rapidly against infection, but is 


durable against severe COVID-19, regardless of variant.*°*° 


Infection with common-cold coronaviruses, and perhaps influenza,** induces only a year-long 
immunity against reinfection, but life-long immunity against severe reinfection.” To what 
extent this pattern reflects waning in biological immunity or immune evasion with virus 
evolution over the global season is unclear. Above results suggest that SARS-CoV-2 may exhibit 
a similar pattern to that of common-cold coronaviruses within few years. Short-term biological 
immunity against reinfection of 3 years may decline as a result of viral evolution and immune 
evasion, leading to periodic (possibly annual) waves of infection. However, the lasting immunity 
against severe reinfection will contribute to a pattern of benign infection. Most primary 
infections would occur in childhood and would likely not be severe. Adults would only 


experience periodic reinfections, also not likely to be severe. 


This study has limitations. We investigated incidence of documented infections, but other 
infections may have occurred and gone undocumented. Undocumented infections confer 


immunity or boost existing immunity, thereby perhaps affecting the estimates (note Section S3). 
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Differences in testing frequency existed between the followed cohorts, but these were small and 
adjusted estimates in sensitivity analyses confirmed similar findings. Gompertz function?’ was 
used to parametrize immunity decay, based on empiric goodness-of-fit, but this analysis serves 
only as an informed exploratory extrapolation that remains to be confirmed with more follow-up 
of cohorts. With Qatar’s young population, our findings may not be generalizable to other 
countries where elderly citizens constitute a larger proportion of the total population. However, 
additional analyses restricting the matched cohorts to those >50 years of age showed findings 


resembling those for the total population. 


Depletion of the primary-infection cohorts by COVID-19 mortality at time of primary infection 
may have biased these cohorts toward healthier individuals with stronger immune responses. 

However, COVID-19 mortality has been low in Qatar’s predominantly young population,!”*° 
totaling 679 COVID-19 deaths (<0.1% of primary infections) up to June 29, 2022, and much 


smaller than the size of the study cohorts. A survival effect seems unlikely to explain or 


appreciably affect study findings, apart perhaps from protection against severe COVID-19. 


As an observational study, investigated cohorts were neither blinded nor randomized, so 
unmeasured or uncontrolled confounding cannot be excluded. While matching was done for sex, 
age, nationality, comorbidity count, and timing of primary infection, this was not possible for 
other factors such as geography or occupation, as such data were unavailable. However, Qatar is 
essentially a city state and infection incidence was broadly distributed across neighborhoods. 
Nearly 90% of Qatar’s population are expatriates from over 150 countries, coming here because 
of employment.'? Most are craft and manual workers working in development projects. !° 


Nationality, age, and sex provide a powerful proxy for socio-economic status in this country.!??* 


27 Nationality alone is strongly associated with occupation. !°?5-27 
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Matching was done to control for factors that affect infection exposure in Qatar.'?*?” The 
matching prescription had already been investigated in previous studies of different 


4,28-31 


epidemiologic designs, and using control groups to test for null effects. These control 


groups included unvaccinated cohorts versus vaccinated cohorts within two weeks of the first 


4,28-30 


dose, when vaccine protection is negligible,*” and mRNA-1273- versus BNT162b2- 


vaccinated cohorts, also in the first two weeks after the first dose.*! These studies have shown 
that this prescription provides adequate control of the differences in infection exposure.*7*3! 


These analyses were implemented using Qatar’s total population with large sample sizes, thus 


minimizing the likelihood of bias. 


In conclusion, protection of natural infection against reinfection wanes and may diminish within 
a few years. Omicron immune evasion accelerates this waning. Meanwhile, protection against 
severe reinfection is very strong with no evidence for waning, regardless of variant, for over 14 


months after the primary infection. 
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Table 1. Baseline characteristics of the eligible and matched primary-infection and infection-naive cohorts in the Pre-Omicron 


Reinfection Study and the Omicron Reinfection Study. 


Pre-Omicron Reinfection Study 


Omicron Reinfection Study 


Full eligible cohorts Matched cohorts” Full eligible cohorts Matched cohorts” 
Characteristics Primary- Infection-naive Primary- Infection- Primary- Infection-naive Primary- Infection- 
infection cohort cohort + | infection cohort naive cohort + infection cohort + infection naive cohort SMD’ 
SMD SMD SMD 
cohort cohort 
N=301,943 N=2,329,039 N=290,638 N=290,638 N=127,075 N=2,329,039 N=120,483 N=120,483 
oe 32 (24-40) 31 (24-39) 0.05! 32 (24-40) 32 (24-40) 0.00! 27 (9-36) 31 (24-39) 0.36! 27 (9-36) 27 (9-36) 0.00! 
Age group 
0-9 years 29,583 (9.8) 258,890 (11.1) 29,022 (10.0) 29,022 (10.0) 32,465 (25.6) 258,890 (11.1) 31,834 (26.4) 31,834 (26.4) 
10-19 years 22,830 (7.6) 168,154 (7.2) 21,900 (7.5) 21,900 (7.5) 11,338 (8.9) 168,154 (7.2) 10,630 (8.8) 10,630 (8.8) 
20-29 years 71,904 (23.8) 585,421 (25.1) 70,011 (24.1) 70,011 (24.1) 27,556 (21.7) 585,421 (25.1) 26,467 (22.0) 26,467 (22.0) 
30-39 years 96,468 (32.0) 736,856 (31.6) 0.07 93,754 (32.3) 93,754 (32.3) 0.00 32,736 (25.8) 736,856 (31.6) 0.40 31,340 (26.0) 31,340 (26.0) 0.00 
40-49 years 52,822 (17.5) 363,820 (15.6) ` 50,356 (17.3) 50,356 (17.3) ` 15,008 (11.8) 363,820 (15.6) i 13,908 (11.5) 13,908 (11.5) : 
50-59 years 20,772 (6.9) 151,131 (6.5) 19,234 (6.6) 19,234 (6.6) 5,453 (4.3) 151,131 (6.5) 4,597 (3.8) 4,597 (3.8) 
60-69 years 5,950 (2.0) 50,033 (2.2) 5,127 (1.8) 5,127 (1.8) 1,937 (1.5) 50,033 (2.2) 1,376 (1.1) 1,376 (1.1) 
70+ years 1,614 (0.5) 14,734 (0.6) 1,234 (0.4) 1,234 (0.4) 582 (0.5) 14,734 (0.6) 331 (0.3) 331 (0.3) 
Sex 
Male 220,978 (73.2) 1,625,533 (69.8) 0.08 212,685 (73.2) 212,685 (73.2) 0.00 85,134 (67.0) 1,625,533 (69.8) 0.06 80,791 (67.1) 80,791 (67.1) 0.00 
Female 80,965 (26.8) 703,506 (30.2) ` 77,953 (26.8) 77,953 (26.8) Í 41,941 (33.0) 703,506 (30.2) i 39,692 (32.9) 39,692 (32.9) i 
Nationality’ 
Bangladeshi 26,984 (8.9) 152,627 (6.6) 25,217 (8.7) 25,217 (8.7) 5,808 (4.6) 152,627 (6.6) 5,280 (4.4) 5,280 (4.4) 
Egyptian 15,517 (5.1) 121,086 (5.2) 15,275 (5.3) 15,275 (5.3) 7,789 (6.1) 121,086 (5.2) 7,525 (6.3) 7,525 (6.3) 
Filipino 23,569 (7.8) 169,647 (7.3) 23,252 (8.0) 23,252 (8.0) 9,244 (7.3) 169,647 (7.3) 9,070 (7.5) 9,070 (7.5) 
Indian 80,738 (26.7) 641,424 (27.5) 80,426 (27.7) 80,426 (27.7) 33,697 (26.5) 641,424 (27.5) 33,122 (27.5) 33,122 (27.5) 
Nepalese 36,149 (12.0) 201,681 (8.7) 33,226 (11.4) 33,226 (11.4) 11,459 (9.0) 201,681 (8.7) 10,373 (8.6) 10,373 (8.6) 
Pakistani 16,779 (5.6) 126,346 (5.4) 0.24 16,182 (5.6) 16,182 (5.6) 0.00 7,591 (6.0) 126,346 (5.4) 0.19 7,228 (6.0) 7,228 (6.0) 0.00 
Qatari 36,177 (12.0) 235,972 (10.1) 36,091 (12.4) 36,091 (12.4) 18,940 (14.9) 235,972 (10.1) 18,813 (15.6) 18,813 (15.6) 
Sri Lankan 9,598 (3.2) 59,805 (2.6) 9,249 (3.2) 9,249 (3.2) 3,190 (2.5) 59,805 (2.6) 2,950 (2.5) 2,950 (2.5) 
Sudanese 8,231 (2.7) 55,207 (2.4) 8,017 (2.8) 8,017 (2.8) 3,638 (2.9) 55,207 (2.4) 3,430 (2.9) 3,430 (2.9) 
Other 
nationalities! 48,201 (16.0) 565,244 (24.3) 43,703 (15.0) 43,703 (15.0) 25,719 (20.2) 565,244 (24.3) 22,692 (18.8) 22,692 (18.8) 
Comorbidity count 
None 241,571 (80.0) 1,991,109 (85.5) 235,153 (80.9) 235,153 (80.9) 104,077 (81.9) 1,991,109 (85.5) 100,533 (83.4) 100,533 (83.4) 
1-2 48,776 (16.2) 268,390 (11.5) 0.15 45,342 (15.6) 45,342 (15.6) 0.00 20,547 (16.2) 268,390 (11.5) 0.15 18,271 (15.2) 18,271 (15.2) 0.00 
3+ 11,596 (3.8) 69,540 (3.0) 10,143 (3.5) 10,143 (3.5) 2,451 (1.9) 69,540 (3.0) 1,679 (1.4) 1,679 (1.4) 


IQR denotes interquartile range and SMD standardized mean difference. 


“Individuals with a documented primary SARS-CoV-2 infection were exact-matched in a 1:1 ratio by sex, 10-year age group, nat 


naive individual. 


ionality, comorbidity count, and calendar week of the SARS-CoV-2 test to the first eligible infection- 


*SMD is the difference in the mean of a covariate between groups divided by the pooled standard deviation. An SMD <0.1 indicates adequate matching. 
#SMD is for the mean difference between groups divided by the pooled standard deviation. 
SNationalities were chosen to represent the most populous groups in Qatar. 
IThese comprise 150 other nationalities in the unmatched primary-infection cohort and 183 other nationalities in the unmatched infection-naive cohort, and 133 other nationalities in the matched primary-infection cohort and 133 
other nationalities in the matched infection-naive cohort in the Pre-Omicron Reinfection Study. These also comprise 145 other nationalities in the unmatched primary-infection cohort and 183 other nationalities in the unmatched 
infection-naive cohort, and 123 other nationalities in the matched primary-infection cohort and 123 other nationalities in the matched infection-naive cohort in the Omicron Reinfection Study. 
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Table 2. Baseline characteristics of the eligible and matched primary-infection and infection-naive cohorts in the COVID-19 


Severity Reinfection Study. 


Full eligible cohorts Matched cohorts” 
Characteristics Primary-infection cohort Infection-naive cohort SMDt Primary-infection cohort Infection-naive cohort SMD* 
N=438,854 N=2,843,496 N=407,214 N=407,214 
Median age (IQR)—years 30 (16-39) 31 (23-39) 0.15? 30 (18-38) 30 (19-38) 0.00# 
Age group 
0-9 years 75,287 (17.2) 324,956 (11.4) 69,059 (17.0) 69,059 (17.0) 
10-19 years 42,913 (9.8) 196,544 (6.9) 36,029 (8.9) 36,029 (8.9) 
20-29 years 97,277 (22.2) 748,686 (26.3) 92,980 (22.8) 92,980 (22.8) 
30-39 years 122,856 (28.0) 880,600 (31.0) 021 117,269 (28.8) 117,269 (28.8) 0.00 
40-49 years 64,930 (14.8) 429,083 (15.1) ` 60,769 (14.9) 60,769 (14.9) f 
50-59 years 25,272 (5.8) 179,111 (6.3) 22,858 (5.6) 22,858 (5.6) 
60-69 years 7,811 (1.8) 64,210 (2.3) 6,492 (1.6) 6,492 (1.6) 
70+ years 2,508 (0.6) 20,306 (0.7) 1,758 (0.4) 1,758 (0.4) 
Sex 
Male 297,317 (67.8) 1,979,916 (69.6) 0.04 280,419 (68.9) 280,419 (68.9) 0.00 
Female 141,537 (32.2) 863,580 (30.4) i 126,795 (31.1) 126,795 (31.1) i 
Nationality’ 
Bangladeshi 29,599 (6.7) 175,765 (6.2) 27,691 (6.8) 27,691 (6.8) 
Egyptian 23,003 (5.2) 136,182 (4.8) 21,534 (5.3) 21,534 (5.3) 
Filipino 32,871 (7.5) 200,002 (7.0) 31,256 (7.7) 31,256 (7.7) 
Indian 108,178 (24.7) 804,623 (28.3) 107,720 (26.5) 107,720 (26.5) 
Nepalese 43,108 (9.8) 253,949 (8.9) 0.26 40,098 (9.9) 40,098 (9.9) 0.00 
Pakistani 22,647 (5.2) 170,868 (6.0) i 21,966 (5.4) 21,966 (5.4) i 
Qatari 71,132 (16.2) 249,740 (8.8) 60,438 (14.8) 60,438 (14.8) 
Sri Lankan 11,912 (2.7) 79,189 (2.8) 11,519 (2.8) 11,519 (2.8) 
Sudanese 11,666 (2.7) 62,039 (2.2) 10,914 (2.7) 10,914 (2.7) 
Other nationalities! 84,738 (19.3) 711,139 (25.0) 74,078 (18.2) 74,078 (18.2) 
Comorbidity count 
None 352,738 (80.4) 2,496,024 (87.8) 334,409 (82.1) 334,409 (82.1) 
1-2 72,153 (16.4) 277,344 (9.8) 0.21 61,720 (15.2) 61,720 (15.2) 0.00 
3+ 13,963 (3.2) 70,128 (2.5) 11,085 (2.7) 11,085 (2.7) 


IQR denotes interquartile range and SMD standardized mean difference. 
“Individuals with a documented primary SARS-CoV-2 infection were exact-matched in a 1:1 ratio by sex, 10-year age group, nationality, comorbidity count, and calendar week of the SARS-CoV-2 test to the first eligible infection- 


naive individual. 


*SMD is the difference in the mean of a covariate between groups divided by the pooled standard deviation. An SMD <0.1 indicates adequate matching. 
#SMD is for the mean difference between groups divided by the pooled standard deviation. 
SNationalities were chosen to represent the most populous groups in Qatar. 


IThese comprise 161 other nationalities in the unmatched primary-infection cohort and 183 other nationalities in the unmatched infection-naive cohort, and 141 other nationalities in the matched primary-infection cohort and 141 


other nationalities in the matched infection-naive cohort. 
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Figure 1. A) Cumulative incidence of infection in the Pre-Omicron Reinfection Study. B) 
Cumulative incidence of infection in the Omicron Reinfection Study. C) Cumulative 
incidence of severe, critical or fatal COVID-19 in the COVID-19 Severity Reinfection 


Study. 
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Table 3. Hazard ratios for the incidence of SARS-CoV-2 infection and incidence of severe, 
critical, or fatal COVID-19 in the Pre-Omicron Reinfection Study, Pre-Omicron 
Reinfection Study, and COVID-19 Severity Reinfection Study. 


Epidemiological measure __Primary-infection cohort — Infection-naive cohort 
Pre-Omicron Reinfection Study H 
Primary outcome 


Total follow-up time (person-weeks) 6,578,466 6,432,430 


Adjusted hazard ratio for SARS-CoV-2 infection’ (95% CT) 0.14 (0.14 to 0.15) 


Effectiveness against severe, critical, or fatal COVID-19* (95% CI) 
Omicron Reinfection Study 


98.0 (95.7 to 99.0) 


Unadjusted hazard ratio for SARS-CoV-2 infection (95% CT) 0.64 (0.62 to 0.66) 


Secondary outcome 


Adjusted hazard ratio for severe, critical, or fatal COVID-19"* (95% CT) 0.11 (0.04 to 0.29) 
i st severe, critical, or fatal COVI 88.6 (70.9 to 95.5) 


Severity Study 


Primary outcome 


Incidence rate of severe, critical, or fatal COVID-19 (per 10,000 person-weeks) 0.01 (0.01 to 0.02) 0.40 (0.36 to 0.44) 


Effectiveness against severe, critical, or fatal COVID-19"' (95% CT) 97.3 (94.9 to 98.6) 


Unadjusted hazard ratio for SARS-CoV-2 infection (95% CI) 0.31 (0.30 to 0.32) 


Effectiveness against SARS-CoV-2 infection in % (95% CD 69.4 (68.6 to 70.3) 
CI denotes confidence interval, COVID-19 coronavirus disease 2019, and SARS-CoV-2 severe acute respiratory syndrome coronavirus 2. 
“Cox regression analysis adjusted for sex, 10-year age group (Table 1), 10 nationality groups (Table 1), comorbidity count (Table 1), and calendar week of the SARS- 
CoV-2 test. 
Severity,” criticality,?? and fatality 


23 were defined according to the World Health Organization guidelines. 
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Figure 2. A) Effectiveness of pre-Omicron primary infection against pre-Omicron 
reinfection. B) Effectiveness of pre-Omicron primary infection against Omicron 
reinfection. C) Effectiveness of primary infection with any variant against severe, critical, 


or fatal COVID-19 due to reinfection with any variant. 
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*Confidence interval could not be estimated using Cox regression because of zero events and was approximated using the confidence interval for the incidence rate ratio. 
tThe negative lower bound for the confidence interval was truncated because the confidence interval was too wide. 
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Figure 3. Extrapolated effectiveness of pre-Omicron primary infection against pre- 


Omicron reinfection, and extrapolated effectiveness of pre-Omicron primary infection 
against Omicron reinfection. 
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Section S1. Laboratory methods and variant ascertainment 


Real-time reverse-transcription polymerase chain reaction testing 


Nasopharyngeal and/or oropharyngeal swabs were collected for polymerase chain reaction 
(PCR) testing and placed in Universal Transport Medium (UTM). Aliquots of UTM were: 1) 
extracted on KingFisher Flex (Thermo Fisher Scientific, USA), MGISP-960 (MGI, China), or 
ExiPrep 96 Lite (Bioneer, South Korea) followed by testing with real-time reverse-transcription 
PCR (RT-qPCR) using TaqPath COVID-19 Combo Kits (Thermo Fisher Scientific, USA) on an 
ABI 7500 FAST (Thermo Fisher Scientific, USA); 2) tested directly on the Cepheid GeneXpert 
system using the Xpert Xpress SARS-CoV-2 (Cepheid, USA); or 3) loaded directly into a Roche 
cobas 6800 system and assayed with the cobas SARS-CoV-2 Test (Roche, Switzerland). The 
first assay targets the viral S, N, and ORF lab gene regions. The second targets the viral N and E- 


gene regions, and the third targets the ORFlab and E-gene regions. 


All PCR testing was conducted at the Hamad Medical Corporation Central Laboratory or Sidra 


Medicine Laboratory, following standardized protocols. 
Rapid antigen testing 


Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antigen tests were performed 
on nasopharyngeal swabs using one of the following lateral flow antigen tests: Panbio COVID- 
19 Ag Rapid Test Device (Abbott, USA); SARS-CoV-2 Rapid Antigen Test (Roche, 
Switzerland); Standard Q COVID-19 Antigen Test (SD Biosensor, Korea); or CareStart COVID- 
19 Antigen Test (Access Bio, USA). All antigen tests were performed point-of-care according to 
each manufacturer’s instructions at public or private hospitals and clinics throughout Qatar with 


prior authorization and training by the Ministry of Public Health (MOPH). Antigen test results 
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were electronically reported to the MOPH in real time using the Antigen Test Management 


System which is integrated with the national Coronavirus Disease 2019 (COVID-19) database. 
Classification of infections by variant type 


Surveillance for SARS-CoV-2 variants in Qatar is based on viral genome sequencing and 
multiplex RT-qPCR variant screening! of random positive clinical samples,”’ complemented by 
deep sequencing of wastewater samples.** Further details on the viral genome sequencing and 
multiplex RT-qPCR variant screening throughout the SARS-CoV-2 waves in Qatar can be found 


in previous publications.” 7-14 
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Section S2. COVID-19 severity, criticality, and fatality classification 


Classification of COVID-19 case severity (acute-care hospitalizations),’° criticality (intensive- 
care-unit hospitalizations),!> and fatality! followed World Health Organization (WHO) 
guidelines. Assessments were made by trained medical personnel independent of study 
investigators and using individual chart reviews, as part of a national protocol applied to every 
hospitalized COVID-19 patient. Each hospitalized COVID-19 patient underwent an infection 
severity assessment every three days until discharge or death. We classified individuals who 
progressed to severe, critical, or fatal COVID-19 between the time of the documented infection 
and the end of the study based on their worst outcome, starting with death,!* followed by critical 


disease,!° and then severe disease. !> 


Severe COVID-19 disease was defined per WHO classification as a SARS-CoV-2 infected 
person with “oxygen saturation of <90% on room air, and/or respiratory rate of >30 
breaths/minute in adults and children >5 years old (or >60 breaths/minute in children <2 months 
old or >50 breaths/minute in children 2-11 months old or >40 breaths/minute in children 1—5 
years old), and/or signs of severe respiratory distress (accessory muscle use and inability to 
complete full sentences, and, in children, very severe chest wall indrawing, grunting, central 
cyanosis, or presence of any other general danger signs)’”.!° Detailed WHO criteria for 


classifying SARS-CoV-2 infection severity can be found in the WHO technical report.!° 


Critical COVID-19 disease was defined per WHO classification as a SARS-CoV-2 infected 
person with “acute respiratory distress syndrome, sepsis, septic shock, or other conditions that 
would normally require the provision of life sustaining therapies such as mechanical ventilation 
(invasive or non-invasive) or vasopressor therapy”.!° Detailed WHO criteria for classifying 


SARS-CoV-2 infection criticality can be found in the WHO technical report. ° 
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COVID-19 death was defined per WHO classification as “a death resulting from a clinically 
compatible illness, in a probable or confirmed COVID-19 case, unless there is a clear alternative 
cause of death that cannot be related to COVID-19 disease (e.g. trauma). There should be no 
period of complete recovery from COVID-19 between illness and death. A death due to COVID- 
19 may not be attributed to another disease (e.g. cancer) and should be counted independently of 
preexisting conditions that are suspected of triggering a severe course of COVID-19”. Detailed 


WHO criteria for classifying COVID-19 death can be found in the WHO technical report.'° 
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Section S3. Additional material for the Results section 
Pre-Omicron Reinfection Study 


The median time of follow-up was 154 days (interquartile range (IQR), 65-224 days) for the 

primary-infection cohort and 151 days (IQR, 61-219 days) for the infection-naive cohort (Figure 
1A). The proportion of individuals who had a SARS-CoV-2 test during follow-up was 29.7% for 
the primary-infection cohort and 36.4% for the infection-naive cohort. The testing frequency was 


0.56 and 0.79 tests per person, respectively. 


The pattern of waning of protection in Figure 2A was fitted to a Gompertz function,” but after 
setting the effectiveness values in months 4-6 after primary infection at 90.5%, the value at the 
7 month. This was done to correct for the likely underestimation of the observed effectiveness 
in these months because of the effect of PCR-positive tests that reflected prolonged infections as 
opposed to true reinfections.!** The fitted Gompertz function suggested that effectiveness 
against reinfection reaches 50% in the 22" month after primary infection, and reaches <10% by 


the 32"! month (Figure 3). 
Omicron Reinfection Study 


The median time of follow-up was 168 days (IQR, 168-168 days) for the primary-infection 
cohort and 168 days (IQR, 147-168 days) for the infection-naive cohort (Figure 1B). The 
proportion of individuals who had a SARS-CoV-2 test during follow-up was 29.9% for the 
primary-infection cohort and 32.5% for the infection-naive cohort. The testing frequency was 


0.50 and 0.55 tests per person, respectively. 


The pattern of waning of protection in Figure 2B was fitted to a Gompertz function,'’ but after 


excluding the effectiveness values before December 1, 2020. Incidence before this date, that is 


medRxiv preprint doi: https://doi.org/10.1101/2022.07.06.22277306; this version posted July 7, 2022. The copyright holder for this preprint 
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 
All rights reserved. No reuse allowed without permission. 


during the first wave, affected mainly the craft and manual workers population of Qatar, that 
constitutes 60% of the total population and where SARS-CoV-2 seroprevalence exceeded 50% 
by the end of the first wave.?*? This population segment has also the lowest testing rates in the 
population.” It is possible that many of these workers may have had mild or asymptomatic 


30.31 that were not documented during the four subsequent waves. !!:!232 This 


reinfections 
undocumented immune boosting, subsequent to the primary infection, may explain the higher- 
than-expected effectiveness for those who had their primary infection prior to December 1, 2020. 


The fitted Gompertz function suggested that effectiveness against reinfection reaches 50% in the 


8" month after primary infection, and reaches <10% by the 15" month (Figure 3). 
COVID-19 Severity Reinfection Study 


The median time of follow-up was 118 days (IQR, 53-235 days) for the primary-infection cohort 
and 111 days (IQR, 52-229 days) for the infection-naive cohort (Figure 1C). The proportion of 
individuals who had a SARS-CoV-2 test during follow-up was 28.0% for the primary-infection 
cohort and 31.7% for the infection-naive cohort. The testing frequency was 0.56 and 0.70 tests 


per person, respectively. 


Since there was no evidence for waning in effectiveness against severe, critical, or fatal COVID- 
19 due to reinfection, there was no relevance to fit a Gompertz function to the effectiveness trend 


after primary infection. 
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Table S1. STROBE checklist for cohort studies. 


pss Recommendation Main Text page 
Title and abstract 1 (a) Indicate the study’s design with a commonly used term in Abstract 
the title or the abstract 
(b) Provide in the abstract an informative and balanced Abstract 
summary of what was done and what was found 
Introduction 
Background/rationale 2 Explain the scientific background and rationale for the Introduction 
investigation being reported 
Objectives 3 State specific objectives, including any prespecified Introduction 
hypotheses 
Methods 
Study design 4 Present key elements of study design early in the paper Methods (‘Study designs and cohorts’) 
Setting 5 Describe the setting, locations, and relevant dates, including Methods (‘Study designs and cohorts’, 
periods of recruitment, exposure, follow-up, and data ‘Cohort matching and follow-up’, ‘Pre- 
collection Omicron Reinfection Study’, ‘Omicron 
Reinfection Study’, & ‘COVID-19 
Severity Reinfection Study’) & Figures 
S1-S3 in Supplementary Appendix 
Participants 6 (a) Give the eligibility criteria, and the sources and methods of | Methods (‘Study designs and cohorts’, 
selection of participants. Describe methods of follow-up “Cohort matching and follow-up’, ‘Pre- 
(b) For matched studies, give matching criteria and number of Omicron Reinfection Study’, ‘Omicron 
exposed and unexposed Reinfection Study’, & ‘COVID-19 
Severity Reinfection Study’) & Figures 
S1-S3 in Supplementary Appendix 
Variables 7 Clearly define all outcomes, exposures, predictors, potential Methods (‘Study designs and cohorts’, 
confounders, and effect modifiers. Give diagnostic criteria, if ‘Cohort matching and follow-up’, ‘Pre- 
applicable Omicron Reinfection Study’, ‘Omicron 
Reinfection Study’, & ‘COVID-19 
Severity Reinfection Study’), Tables 1- 
2, & Sections S1& S2 in 
Supplementary Appendix 
Data sources/ 8* For each variable of interest, give sources of data and details of | Methods (‘Study population and data 
measurement methods of assessment (measurement). Describe comparability | sources’ & ‘Statistical analysis’), 
of assessment methods if there is more than one group Tables 1-2, & Sections S1 & S2 in 
Supplementary Appendix 
Bias 9 Describe any efforts to address potential sources of bias Methods (‘Cohort matching and follow- 
up’) 
Study size 10 Explain how the study size was arrived at Figures S1-S3 in Supplementary 
Appendix 
Quantitative 11 Explain how quantitative variables were handled in the Methods (‘Cohort matching and follow- 
variables analyses. If applicable, describe which groupings were chosen up’) & Tables 1-2 
and why 
Statistical methods 12 (a) Describe all statistical methods, including those used to Methods (‘Statistical analysis’) 
control for confounding 
(b) Describe any methods used to examine subgroups and Methods (‘Statistical analysis’, 
interactions paragraph 3) 
(c) Explain how missing data were addressed NA, see Methods (‘Study population 
and data sources’) 
(d) If applicable, explain how loss to follow-up was addressed NA, see Methods (‘Study designs and 
cohorts’, paragraph 1) 
(e) Describe any sensitivity analyses Methods (‘Statistical analysis’, 
paragraph 3) 
Results 
Participants 13* (a) Report numbers of individuals at each stage of study—eg Figures S1-S3 in Supplementary 
numbers potentially eligible, examined for eligibility, Appendix 
confirmed eligible, included in the study, completing follow- 
up, and analysed 
(b) Give reasons for non-participation at each stage 
(c) Consider use of a flow diagram 
Descriptive data 14 (a) Give characteristics of study participants (eg demographic, Tables 1-2 
clinical, social) and information on exposures and potential 
confounders 
(b) Indicate number of participants with missing data for each NA, see Methods (‘Study population 
variable of interest and data sources’) 
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(c) Summarise follow-up time (eg, average and total amount) Results (‘Pre-Omicron Reinfection 
Study’, ‘Omicron Reinfection Study’, 
& ‘COVID-19 Severity Reinfection 
Study’), Figure 1, & Table 3 
Outcome data 15 Report numbers of outcome events or summary measures over | Results (‘Pre-Omicron Reinfection 
time Study’, ‘Omicron Reinfection Study’, 
& ‘COVID-19 Severity Reinfection 
Study’), Figure 1, & Table 3 
Main results 16 (a) Give unadjusted estimates and, if applicable, confounder- Results (‘Pre-Omicron Reinfection 
adjusted estimates and their precision (eg, 95% confidence Study’, ‘Omicron Reinfection Study’, 
interval). Make clear which confounders were adjusted for and | & ‘COVID-19 Severity Reinfection 
why they were included Study’), & Table 3 
(b) Report category boundaries when continuous variables Tables 1-2 
were categorized 
(c) If relevant, consider translating estimates of relative risk NA 
into absolute risk for a meaningful time period 
Other analyses 17 Report other analyses done—eg analyses of subgroups and Results (‘Pre-Omicron Reinfection 
interactions, and sensitivity analyses Study’, ‘Omicron Reinfection Study’, 
& ‘COVID-19 Severity Reinfection 
Study’) & Figures 2-3 
Discussion 
Key results 18 Summarise key results with reference to study objectives Discussion, paragraphs 1-4 
Limitations 19 Discuss limitations of the study, taking into account sources of | Discussion, paragraphs 5-8 
potential bias or imprecision. Discuss both direction and 
magnitude of any potential bias 
Interpretation 20 Give a cautious overall interpretation of results considering Discussion, paragraph 9 
objectives, limitations, multiplicity of analyses, results from 
similar studies, and other relevant evidence 
Generalisability 21 Discuss the generalisability (external validity) of the study Discussion, paragraphs 5-7 and Table 
results S2 in Supplementary Appendix 
Other information 
Funding 22 Give the source of funding and the role of the funders for the Sources of support & 
present study and, if applicable, for the original study on which | acknowledgements 
the present article is based 


Figure S1. Cohort selection for investigating immune protection of primary infection with a pre-Omicron variant against re- 


infection with a pre-Omicron variant (Pre-Omicron Reinfection Study). 


419,963 Individuals with a documented primary SARS-CoV-2 
infection between February 28, 2020 and November 30, 2021 


118,020 Excluded 

e 33,686 Had a vaccination record before the 
documented primary infection 

e 70,993 Had a vaccination record before the 
start of the follow-up (90 days after the 
documented primary infection) 

e 12,603 Had <90 days after the documented 
primary infection 

e 738 Died before the start of the follow-up, of 
which 523 were COVID-19-related deaths 


301,943 Individuals with a documented 
primary infection and no vaccination record 


11,305 Not matched 


290,638 Individuals with a documented primary infection were 
matched” to individuals who were alive, infection free, and 
unvaccinated as of the start of follow-up (90 days after the 

documented primary infection) 


v 


1,806 Reinfected with SARS-CoV-2 
o 433 PCR-confirmed symptomatic infections* 
= 5 progressed to severe COVID-19 
o 1,373 PCR-confirmed non-symptomatic infections* 
= | progressed to severe COVID-19 
=" 1 progressed to fatal COVID-19 
e 110,266 Censored because they received Dose | 


e 104,310 Censored because their match received Dose | 
83 Deaths during follow-up, of which 13 were COVID-19- 
related deaths related to the primary infection and not due to a 
reinfection= 
74,173 Followed until end of study (November 30, 2021) 


2,927,678 Individuals with PCR-negative SARS-CoV-2 
tests between February 5, 2020 and November 30, 2021 


598,639 Excluded 
e 432,191 Had a vaccination record before 
the PCR-negative test 
© 166,353 Had a record of prior infection 
before the PCR-negative test 
e 95 Had no ascertained death date 


2,329,039 Uninfected individuals with no 
vaccination record 


2,038,401 Not matched 


290,638 Uninfected control individuals with no vaccination 
record were matched" to individuals with a documented primary 
infection. Controls had to be alive, infection free, and 
unvaccinated as of the start of follow-up (90 days after the 
documented primary infection of their match) 


11,957 Infected with SARS-CoV-2 
o 4,011 PCR-confirmed symptomatic infectionst 
= 199 progressed to severe COVID-19 
= 10 progressed to critical COVID-19 
= 9 progressed to fatal COVID-19 
o 7,946 PCR-confirmed non-symptomatic infections‘ 
= 98 progressed to severe COVID-19 
= 9 progressed to critical COVID-19 
= 3 progressed to fatal COVID-19 
99,966 Censored because they received Dose | 
106,002 Censored because their match received Dose 1 
175 Non-COVID-19-related Deaths during follow-up 
72,538 Followed until end of study (November 30, 2021) 


COVID-19 denotes coronavirus disease 2019, PCR polymerase chain reaction, and SARS-Co¥V-2 severe acute respiratory syndrome coronavirus 2. 

*Cohorts were exact-matched in a 1:1 ratio by sex, 10-year age group, nationality, comorbidity count, and calendar week of the SARS-CoV-2 test. 

tA PCR-confirmed symptomatic infection was defined as a PCR-positive test conducted because of clinical suspicion due to presence of symptoms compatible with a respiratory tract infection. A PCR-confirmed non-symptomatic infection was defined 
as an infection diagnosed for a reason for PCR testing other than clinical suspicion. Some of the latter infections developed symptoms after diagnosis and some progressed to severe forms of COVID-19. 

These deaths occurred among persons who had the death >90 days after primary infection but had no record of a documented reinfection during follow-up. 
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Figure S2. Cohort selection for investigating immune protection of primary infection with a pre-Omicron variant against re- 


infection with Omicron subvariants (Omicron Reinfection Study). 


419,963 Individuals with a documented primary SARS-CoV-2 2,927,678 Individuals with negative SARS-CoV-2 tests 
infection between February 28, 2020 and November 30, 2021 between February 5, 2020 and November 30, 2021 


292,888 Excluded 


¢ 33,686 Had a vaccination record before the 598,639 Excluded 

documented primary infection e 432,191 Had a vaccination record before the 
e 257,517 Had a vaccination record before the SARS-CoV-2-negative test 

start of the follow-up" 166,353 Had a record of prior infection before 
¢ 898 Died before the start of the follow-up, of the SARS-CoV-2-negalive test 

which 550 were COVID-19-related deaths o 95 Had no ascertained death date 


è 787 Had a record of reinfection before the 
start of the follow-up 


127,075 Individuals with a documented 2,329,039 Uninfected individuals with no 
primary infection and no vaccination record vaccination record 


6,592 Not matched 2,208,556 Not matched 


120,483 Individuals with a documented primary infection were 120,483 Uninfected control individuals with no vaccination 
matched! to individuals who were alive, infection free, and record were matched’ to individuals with a documented primary 
unvaccinated as of the start of follow-up” infection. Controls had to be alive, infection free, and 
unvaccinated as of the start of follow-up” 


7,995 Reinfected with SARS-CoV-2 
o 1,257 PCR-confirmed symptomatic infections* 
= 3 progressed to severe COVID-19 
o 3,270 PCR-confirmed non-symptomatic infections* 
= | progressed to severe COVID-19 
o 3,468 Rapid-antigen-confirmed infections 
= | progressed to severe COVID-19 
7,630 Censored because they received Dose 1 
6,026 Censored because their match received Dose | 
27 Non-COVID-19-related deaths during follow-up 
98,805 Followed until end of study (June 5, 2022) 


è 12,230 Infected with SARS-CoV-2 
o 2,096 PCR-confirmed symptomatic infections* 
= 8 progressed to severe COVID-19 
= 6 progressed to critical COVID-19 
= 3 progressed to fatal COVID-19 
o 4,954 PCR-confirmed non-symptomatic infections* 
= 10 progressed to severe COVID-19 
= | progressed to critical COVID-19 
o 5,180 Rapid-antigen-confirmed infections 
= 8 progressed to severe COVID-19 
= 2 progressed to fatal COVID-19 
5,881 Censored because they received Dose | 
7,376 Censored because their match received Dose 1 
33 Non-COVID-19-related deaths during follow-up 
94,963 Followed until end of study (June 5, 2022) 


COVID-19 denotes coronavirus disease 2019, PCR polymerase chain reaction, and SARS-CoV-2 severe acute respiratory syndrome coronavirus 2. 

“Follow-up for cases was from 90 days after the documented primary infection if start of follow-up coincided with the SARS-CoV-2 Omicron wave which started on December 19, 2021, or from December 19, 2021 for earlier infections. Follow-up for controls was 
from 90 days after the documented primary infection of their match if start of follow-up coincided with the SARS-CoV-2 Omicron wave which started on December 19, 2021, or from December 19, 2021 for earlier infections 

*Cohorts were exaci-matched in a 1:1 ratio by sex, 10-year age group, nationality, comorbidity count, and calendar week of the SARS-CoV-2 test. 

+A PCR-conlirmed symptomatic infection was defined as a PCR-positive test conducted because of clinical suspicion due to presence of symptoms compatible with a respiratory tract infection. A PCR-conlirmed non-symptomatic infection was defined as an 
infection diagnosed for a reason for PCR testing other than clinical suspicion. Some of the latter infections developed symptoms after diagnosis and some progressed to severe forms of COVID-19. There was no record for reason for testing and presence of 
symptoms for most of the rapid antigen testing and thus we did not report the symptom status for infections diagnosed using rapid antigen testing, 
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Figure S3. Cohort selection for investigating immune protection of primary infection with any variant against severe, critical, 
or fatal COVID-19 due to reinfection with any variant. 


808,396 Individuals with a documented primary SARS-CoV-2 3,647,782 Individuals with negative SARS-CoV-2 tests 
infection between February 28, 2020 and June 5, 2022 between February 5, 2020 and June 5, 2022 


369,542 Excluded 
e 289,329 Had a vaccination record before the 804,286 Excluded 
documented primary infection e 576,359 Had a vaccination record before 
e 75,362 Had a vaccination record before the the negative test 
start of the follow-up e 227,827 Had a record of prior infection 
e 4,009 Had <90 days after the documented before the negative test 
primary infection e 100 Had no ascertained death date 
e 842 Died before the start of the follow-up, of 
which 565 were COVID-19-related deaths 


438,854 Individuals with a documented 2,843,496 Uninfected individuals with no 
primary infection and no vaccination record vaccination record 
—>| 2,436,282 Not matched 


407,214 Individuals with a documented primary infection were matched’ to 407,214 Uninfected control individuals with no vaccination record were 

individuals who were alive, infection free, and unvaccinated as of the start of matched’ to individuals with a documented primary infection. Controls had to 

follow-up (90 days after the documented primary infection) be alive, infection free, and unvaccinated as of the start of follow-up for their 
match (90 days after the documented primary infection of their match) 


31,640 Not matched 


7,082 Reinfected with SARS-CoV-2 21,645 Infected with SARS-CoV-2 
o 9 progressed to severe COVID-19 o 315 progressed to severe COVID-19 
o | progressed to fatal COVID-19 o 25 progressed to critical COVID-19 


118,301 Censored because they received Dose 1 o 18 progressed to fatal COVID-19 
111,423 Censored because their match received Dose | 107,062 Censored because they received Dose 1 
112 Deaths during follow-up, of which 12 were COVID-19-related 113,903 Censored because their match received Dose | 
deaths due to primary infection and not due to a reinfection’ 184 Non-COVID-19-related deaths during follow-up 
e 170,296 Followed until end of study (June 5, 2022) 164,420 Followed until end of study (June 5, 2022) 


COVID-19 denotes coronavirus disease 2019 and SARS-CoV-2 severe acute respiratory syndrome coronavirus 2. 
‘Cohorts were exact-matched in a 1:1 ratio by sex, 10-year age group, nationality, comorbidity count, and calendar week of the SARS-CoV-2 test. 
*These deaths occurred among persons who had the death >90 days after primary infection but had no record of a documented reinfection during follow-up. 
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Table S2. Representativeness of study participants. 


a Category e er a OTTO nT nT teCyE eee 
Disease, problem, or condition under Protection conferred by natural SARS-CoV-2 infection against SARS-CoV-2 
7 inyestigation leee ENfEction and against severe, critical, or fatal COVID-19 due to reinfection. 
= Special considerations related tO nnnnnnnininnnnnnnntnnnninnninnnnnnnnnnnn 
Sex and gender Three national, matched, retrospective, target-trial cohort studies were conducted 


to compare incidence of SARS-CoV-2 infection and COVID-19 severity among 
those unvaccinated but with a documented SARS-CoV-2 primary infection, to 
incidence among those infection-naive and unvaccinated. Cohorts were exact- 
matched by sex to control for potential differences in the risk of exposure to 
ee E SARS-CoV-2 infection by SeX. l l 
Age Cohorts were exact-matched by 10-year age group to control for potential 
differences in the risk of exposure to SARS-CoV-2 infection by age. Nonetheless, 
with the young population of Qatar, our findings may not be generalizable to other 
countries where elderly citizens constitute a larger proportion of the total 
population. Having said so, the additional analyses restricting the matched cohorts 
Race or ethnicity group Cohorts were exact-matched by nationality to control for potential differences in 
the risk of exposure to SARS-CoV-2 infection by nationality. Nationality is 
oe associated with race and ethnicity in the population of Qatar. le 
Geography Individual-level data on geography were not available, but Qatar is essentially a 
city state and infection incidence was broadly distributed across the country’s 
neighborhoods/areas. Cohorts were exact-matched by nationality to control for 
potential differences in the risk of exposure to SARS-CoV-2 infection by 
nationality. Qatar has unusually diverse demographics in that 89% of the 
population are international expatriate residents coming from over 150 countries 


may have (partially) controlled the differences in occupational risk, in 
consideration of the association between nationality and occupation in Qatar. To 
ensure that all individuals in all cohorts have a record of active residence in Qatar 
at the same calendar time, individuals who were diagnosed with a primary 
infection in a specific week were matched to infection-naive individuals who had a 
cine TECO of a SARS-CoV-2 negative test in that same calendar week. llle 

Overall representativeness of this study | The study was based on the total population of Qatar and thus the study population 
is broadly representative of the diverse, by national background, but young and 
predominantly male, total population of Qatar. While there could be differences in 
the risk of exposure to SARS-CoV-2 infection by sex, age, nationality, and 
comorbidity, cohorts were exact-matched by these factors to control for their 
potential impact on our estimates. Given that only 9% of the population of Qatar 
are >50 years of age, our estimates may not be generalizable to other countries 
where elderly citizens constitute a larger proportion of the total population. 
COVID-19 coronavirus disease 2019, and SARS-CoV-2 severe acute respiratory syndrome coronavirus 2. 
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